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Abstract 

The technique of muon mean energy determination in multilayer detectors 
is developed. The mean energy is measured by means of average small bursts m 
i.e. the number of electrons and positrons generated by muons in the detecting 
layers of device via three basic processes — creation of e + e~ pairs, 5-electrons 
and bremsestrahlung. The accuracy of the method is considered. 
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1 Introduction. 

For investigation of the penetrating component of cosmic radiation by 
means of multilayer detectors it is essential usually to estimate the energy of 
registered events. This task in cosmic ray physics faces the serious technical 
obstacles, especially for large detectors and for muons of very high energies. 
However many problems in cosmic ray physics can be solved if one restricts to 
the estimation of the mean energy of muons in samples of events of different 
types, e.g. in muon groups of some fixed multiplicity n. The method of 
muon mean energy estimation is evaluated below. 

This method is based on the analysis of a small bursts accompanying the 
passage of penetrating particles through the layers of filter and detecting 
elements of device. These bursts are generated by muons in the filter layers of 
device via three basic processes — creation of e + e~ pairs, (5-electron emission 
and bremsestrahlung. The value of mean burst m turned out to be very 
effective parameter to evaluate the mean energy of penetrating particles in 
multilayer detectors. 

The high efficiency of penetrating particle registration and energy esti- 
mation is due to the circumstances quoted below. Despite the fact that the 



muon energy losses caused by the pair production process are close to those 
caused by the bremsestrahlung, the cross-section of the first-mentioned pro- 
cess is essentially higher. It increases with the muon energy value and at 
E > 10 12 eV the probability of e + e" pair creation in the absorber layer of the 
optimal thickness is close to the unit. The thickness of the absorber layers 
is to be chosen with the account of the requirement of the minimization of 
particle number fluctuations in avalanches. The calculations [[If show that 
at the muon energies of the order of 10 12 -r- 10 13 eV the minimal dispersion is 
ensured with the thickness of lead absorber of 10 t-units, i.e. ~ 5cm. 

The approach stated below becomes most effective for the investigation of 
high energy muons by means of visual track method of observation. Visual 
observation allows to register even the smallest (single particle) bursts, to 
distinguish easily the events by their nature and to handle the rare events of 
high multiplicity. The influence of transition effects are simple to be taken 
into account since it is easy to separate the particles created in absorber from 
those created in the lid and the walls of spark chambers. The application of 
the multilayer system of spark chambers to the evaluation of muon energy 
was proposed earlier in j2j]. 

Up to day in our experiments (e.g. j3j]) we used the version of this ap- 
proach described in [f3j. In the present work some further development and 
verification of the method is proposed. 

For the analysis of the approach the really working |3j multilayer spark 
detector is considered below. The detector consists of 8 layers of spark 
chambers separated by a 5cm thick lead absorbers. The area of each layer 
of the detector is 4.8m 2 . The spark chambers are made of glass and contain 
neon of high purity at atmospheric pressure. The electrodes are made of 
1mm duraluminium. The detector is located under the rock at the depth 
of 190 hg ■ cm~ 2 . This determines the threshold energy of registered muons 
- about 35GeV. The method proposed is in principle applicable for the 
systems containing the registering elements of any type (e.g. Heiger-Miiller 
counters, neon tubes etc.) 

In chapter 2 the main parameters used are considered. In chapters 3 and 
4 the principle of muon mean energy estimation in multilayer detectors is 
justificated. Chapter 5 is dedicated to the questions of device calibration. 
In chapter 6 actually the process of muon energy estimation is described. 
The dependence E(m) is given for the detector under consideration. 



2 



2 The main input parameters. 



In the real observations of high and superhigh energy cosmic ray muons the 
following parameters are usually available or easily measurable: 

1. The mean threshold energy T = T{x) of registered muons. This value 
depends on the thickness x of the filter located over detector (rock, 
water etc.) and its absorbing features. 

2. The exponent (3 of the spectrum of accompanying bursts Y{m) oc 
m -(/ 3 + 1 ) - According to |5j], in the case of power form of muon energy 
spectrum P(E) oc E~^ 1+1 ^ and with the condition of equilibrium of ac- 
companying electromagnetic cascade ( i.e. in the vicinity of avalanche 
maximum) the exponents (3 and 7 approximately coincide. 

3. The value of mean burst m (i.e. the mean number of e + and e~ ac- 
companying the passage of penetrating particle through the multilayer 
system per one layer). 

Note that the value of m in the most of cases is determinable with much 
less error than T and (3. 

If the differential energetic spectrum of penetrating particles registered 
under the filter is defined in normalized power form 

7 / e r l7+1) 

then their mean energy (under the filter) 

00 

<E> = jE-P{E)dE (2.2) 


is equal to 

<E> = 7 >1. (2.3) 

7-1 



Now, if in fl2T3Q one accepts for valuations 7 = (3 (the conditions of such 
substitution correctness were given above in item 2), it becomes possible to 
determine the mean energy valuation E through the energetic threshold T 
and the exponent (3 of the differential spectrum of small bursts. However ac- 
tually the conditions necessary for such application of ( |2.3|) are not fulfilled. 



Even for the single muons the spectrum differs from (|2.1| ); it is significantly 



3 



flatter at low energies. On the other hand, for some samples of events it 
turns out that (3 < 1. In this case the mean energy < E > formally is 
not defined at all. However m in these samples is steadily recurring in the 
consecutive experiments, so in any case one can consider the value of mean 
burst rn as the well defined one. This means that the burst spectrum is also 
not of pure power-type. It gets steeper at the big bursts. 



3 Determination of mean energy of muons. 

The parametrization of muon energy (Eq) spectra in atmosphere (both for 
single muons and for those which are a component of the events of any other 
type) will be taken in the normalized form 



P{Eo) -%- 2Fl (i,i;i + i;i-T /K ) ■ (31) 

Here 2^1 is the Gaussian hypergeometric function, Tq - - the characteristic 
energy describing flattening of the spectrum at low energies (of the order of 
lGeV); Kq - - the critical energy characterizing the competition of two pro- 
cesses — decay of pions and their interaction with air nuclei at the height of 
effective generation of muons j7j] . Generally speaking, the analogous process 
for i^-mesons is to be taken into account as well. For this reason in ( |3.1|) 
it is taken Kq ~ lOOGeV - - the effective critical energy approximately de- 
scribing the evolution of the flux of tt- and i^-meson mixture |jj . For a not 
very large values of filter thickness x it is possible not to take into account 
the large fluctuations of muon energy losses and thus the mean energy losses 
may be taken in the standard form 

dE 

- — = a + bE, (3.2) 

dx 

where a and b are the parameters characterizing the material of the filter. 
In this case the form of the spectrum is conserved through the evolution of 
spectrum down the filter depth |8j , while the initial energy Eq of muons (the 
energy at the entry to the filter) is the integral of evolution. It is connected 
with the muon energy E under the filter by the relation 

E = (C + E)e bx -C, (3.3) 

where C = a/b is the parameter (energy dimensioned) characterizing the 
material of filter (for the small depths x and small energies E one obtains 
the usually used approximate expression Eq ~ E + ax ) . 
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The relation ( |3T3|) describing the full evolution effect of the spectrum with 
respect to the depth x makes it possible to re-write the spectrum of muons 
under the filter in the form dependent on E and analogous to 



P(Elx) = _JL_ . (1 + TOP ■ (1 + E/Kjx))- 1 
(l > T(x) 2^(1,1; 7 +l;l-T(x)/if(.T)) ' 1 ' ' 

Here the functions T{x) and K(x) describe the evolution of the characteristic 
energy Tq and critical energy Kq with the depth x of the filter: 

T T(x) = C-(C- T )e- bx ~ T + (a - T b)x; 1 
K = K(x) = C-(C- K Q )e- bx ~K + (a- K Q b)x. J [ ' } 

Further we will accept that the energy spectra in the samples of different 
kinds of events (in particular, differing by the number of muons in muon 
groups) may be described by means of functions of ( p. 4|) -class. Here we 
take the natural assumption that energetic parameters T and K specified 
by the process of absorption in filter and by the n- and i^-decays at high 
altitudes are the same for all types of samples of events. We assume that 
the difference of spectra of these samples is described completely by the 
difference of values of their exponents 7. Now, once evaluating the value of 
7 by means of any method one can estimate the mean energy of muon by 
the distribution 



()> 7-1 2 F 1 (i ) i; 7 +i;Z)' ( ' 

where 

Z = Z{x) = 1 - T(x)/K{x). (3.7) 

4 Evaluation of the effective exponent of the spec- 
trum. 

To evaluate the effective value of 7 we use the value of mean burst m which is 
easy to obtain for the given sample of events. On the basis of [jl], £j the mean 
electromagnetic cascade accompanying the passage of muon of fixed energy 
E through the layer of the absorber of the optimal thickness was numerically 
calculated. The obtained E'-dependence of m can be approximated as 

fh{E) = B-{E/e) a , (4.1) 
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where ffi is an average burst at fixed energy E of the muon; e = lGeV - - is 
the energy scale; the parameters a, B are determined by the properties and 
the size of the absorber in the detector. The exponent a one can estimate 
through the numerical calculations with sufficiently high accuracy. The value 
of the normalization parameter B one has to estimate through calibration 
of the device by the well known flux of a single cosmic ray muons (see 
further). Of course, the more accurate results can be achieved by measuring 
the relation like ( |4.1|) with a monochromatic beam of muons from accelerator 
at a different energies. However, only the first way is practically available 
in cosmic ray experiments. 

Averaging the approximate relation ( |4.1| ) for the m value over the spec- 
trum P(E\x) of muon energies (|5.4|) one obtains for the twice-averaged burst 
to be used further again: 




m =< m(E) > = J rn{E') ■ P{E'\x)dE' = 

o 

r(a + i)r(7-a) 2 *i(i,a + i;7 + i;£) 



r(7) 2 F 1 (1,1; 7 +1;Z) 
or, writing this relation in different way 



(4.2) 



$(7, m\a, B, T, K) - 
B (T\ a r(a+l)r(7-a) 2^1(1, a + 1; 7 + 1; Z) 



m \ej r( 7 ) 2 F 1 (1,1; 7 +1;Z) 



1 = 0. (4.3) 



Equation (|4.3| ) defines the implicit function j(m\a, B 1 T, K). This equa- 
tion is to be solved numerically preliminarily being brought to the form 
convenient for the iteration method: 

1 + G . r(7-°) . »fitt,°+i;7 + i;Z) x (44) 

7 r( 7 -l) 2 F 1 (l,l; 7 +l;Z) ' 7 ( ' 

Here the value 

73 / r T\ a 

G = -■(-) -r(a + l) (4.5) 
m \e J 

contains the full experimental information on the sample under investiga- 
tion. The properties of the detector are reflected both in value of G and in 
the parameters a and Z (see Q377D ). The iteration process fl§21 ) converges 
rapidly. It is convenient to start with the value 7 = 2.0. 
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The exponent 7 evaluated from (|4.4|) being inserted in ( p.fc>|) defines the 



searched quantity of the single muon mean energy in the sample under in- 
vestigation. 

So the task is naturally divided into two stages. The first one is calibra- 
tion of the device, i.e. the procedure of evaluation of the basic parameters 
(a, £>, T andK). The second step is application of calibrated device actually 
to the main task — estimation of separate muon mean energy in the sample 
of penetrating events. 

5 Calibration of device. 

As one can easily see from above, the goal of our calculations - - a pair of 
quantities 7 and E =< E > - - is completely fixed by the measurable quan- 
tity m =< fh > and by the set of device parameters represented below by 
the algebraic vector / = (a, B, T, K). Further the description of the method 
will be accompanied by numerical illustration applied to our installation (see 
"Introduction"). 

Parameter a in (|4. 1| ) is obtained through calculations and does not de- 
pend on other parameters. Parameters T and K are obtained from 



They do correlate strongly and are defined by properties of the filter above 
the installation and the values of Tq and Kq quantities for atmospheric 
muons. 

The parameter of proportionality B is estimated from the precise value 
of measured mean burst m\ in events with single cosmic ray muons. The 
spectrum C |5TT| ) of single muons in atmosphere is well known and has the 
exponent value 71 = 2.65 ± 0.05. The calibration measurements of mean 
accompanying burst valuation rh\ for the threshold energy T = 35 ± 2 GeV 
of registered single muons have been carried out (the filter thickness is known 
with accuracy of ~ 5%). For the sample of 150 000 single cosmic ray muons 
the value of mean burst for the considered detector appeared to be rh\ = 
0.173 ± 0.002 (the error is statistical only). 

From fljOj ) it follows: 



B - 5(mi,7i,a,T, K) = 

m ( T Y" r ^'> 2^(1,1,71 + 1,2) (51) 

l \e) r(a + l)r( 7 i - a) 2 F 1 (l,a + l, 7 i + l,Z\ 1 ' ' 

Thus the vector of installation parameters f(a,B,T, K) is to be calculated 
through the vector of initial parameters 
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p = (b,C,x,T , X ,a,7i,mi), (5.2) 

which were defined above. 

Since all the eight parameters of p vector are to be defined from indepen- 
dent experiments it seems to be reasonable to accept the approach to these 
valuations as to normally distributed independent random quantities. For 
the particular case of our installation we accept: 





' Pi 


b = 


(4.0 ± 0.4) • 10" 6 cm 2 • g- 1 


[10] 




P2 


C = 


460.0 ±46.2 GeV 


m 




P3 


x = 


19000.0 ±950.0 g-cm- 2 


H 


p = < 


Pa 


f = 


1.50 ±0.15 GeV 


H (5.3) 
i 


P5 


Ko = 


110.0 ±5.50 GeV 




P6 


a = 


0.90 ±0.05 






P7 


7i = 


2.65 ±0.05 


i 




. P8 


rhi = 


0.173 ±0.002 




The transition fro 


m initial parameters p to consolidated parameters f(p) 



for the normally distributed quantities p with covariation matrix valuation 
V (in this case diagonal one) is reduced to calculation of valuation 



IJ 



/i = l, 2,3,4. 



(5.4) 



dpidpj 

Here a correction proportional to the second derivatives ([5T4]) of f(p) is small, 
so all such corrections are further neglected. The covariation matrix of the 
valuations of consolidated parameters 



D 



jJLV 



= E v l3 



dfu 



dpi 



dpj 



(5.5) 



naturally is nondiagonal one. Here at least some of the derivatives above are 
to be calculated numerically. The quantities of valuations of consolidated 
parameters and of matrix D are given in Table 1. 

Besides, in Fig.l the dependence B(a) is given for the device under con- 
sideration. Here the error for a is taken to be zero. The dashed lines are 
the borders of one standard deviation. 

According to ( |3.fcjQ one obtains for cosmic ray single muons at the thresh- 
old energy f = 35 GeV the value E x = (29.8±1.9) GeV. Note that valuation 
of E at T — 35 GeV directly from the experimental data of highly precise 
work [[TTj leads to the value E[ = (32.9 ± 2.0) GeV. 
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6 Evaluation of the exponent 7 and of the mean en- 
ergy E of muons from the sample. 



Let us accept that for a sample of penetrating events the experimentally 
determined mean burst m is a normally distributed random quantity. This 
is approximately valid for sufficiently large samples. In such case when one 
obtains the valuations of 7 and E through ( p~3[) and (|3T^|) . It is possible 
to consider them as a substitution of random variables analogous to ( |b.4|) 
and (|5T51) . Neglecting the corrections to mean values (as we did above) at 
calculation of 7 we restrict ourself to the solution of equation (fO|) (the 
error of calculation must be smaller than statistical one for 7 — see below). 
Namely this value of 7 we will use to obtain the valuation of E through the 
relation ( |3.tj|) . Since the vector of consolidated parameters / and valuation 
of mean burst m are statistically mutually independent, it is easy to obtain 
the following formulae for the dispersion of 7 and for covariation of 7 with 



4 

A( 7 ) = £ *V ■ — 



df, 



+ <7, 



f,rh 



4 . 

A(7,/p)= E D^- —j- 



i2 



J f,m 



dm 
.dfp 

f,m f,m 

p = l,2,3,4. 



(6.1) 



Here <7^ is the valuation of dispersion of the mean burst m which is got 
from the experiment. Derivatives of 7(7, n) are calculated from (|4.5|) as a 
derivatives of implicit function: 



d7 d$/df v 



$7 
dm 



v 



d^/dm 
1,2,3,4. 



(6.2) 



The dispersion of mean energy valuation for muons of the sample is obtained 
analogously {E = 7), see (|3~6|) ) : 



\l,V=\ U J jJL 



dE 



df, 



+ 



/,7 
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4 ~ dE 
+2-EA( 7 ,/; 



+ 

/,7 



+A( 7 ) 



<9£ 

S7 



2 



(6.3) 



For example the sample of 4-muon events with the mean burst rh^ = 0.45 ± 
0.07 leads to the valuations of 74 and E4. 

74 = 1.65 ±0.13 , E 4 = 91 ± 21GeV. (6.4) 

In Fig. 2 the set of dependencies E{m) is given for the set of values a = 
0.80; 0.85; 0.90; 1.00 — this parameter may vary most probably in different 
devices (if one pais no attention to the possibility of soil-thickness variation 
above the installation). Here a m = 0, a a = is accepted. It is easy to see 
that in the region of relatively small energies (bursts) the result is not very 
sensitive to the variations of a. 

In Fig. 3 the dependencies 7(771) and E(m) are given for the above quoted 
values of consolidated parameters / and for a m = 0, i.e. the corridor of stan- 
dard deviations indicated for central values of mean bursts is the systematic 
error of the detector. Separately standing points on Fig. 3 are the following: 

1. - - direct evaluation of exponent $4 in sample of 4-muon events for the 
case of exponential cutoff of the spectrum in the area of large bursts. 



value E4, calculated according to ( |3.6| ) for the above-mentioned value 



of /?4 by identification 74 = 04 (see "Introduction"). 

- direct evaluation of exponent in sample of 4-muon events for 
purely power spectrum of bursts (argument m is calculated from this 
spectrum). 

- value £4 corresponding to the equality 74 = j3' A is calculated accord- 
ing to fl2T3p . 



Only statistical errors are given for these four points, as well as for the 
central point, got by the regular method proposed. 

Note that evaluation of E is rather stable in respect to variation of spec- 
trum shape at large bursts. In particular if one takes for spectrum in the 
sample the power approximation (|2.1| ) but estimates m neglecting the very 
large bursts, the values of E obtained for close ra-s differ insignificantly. 
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7 Conclusion. 



The proposed method for muon mean energy estimation with the help of 
accompanying burst mean value m of the sample of penetrating events is 
based on the solution of implicit equation (|0|) for the spectrum exponent 
7. The shape of this equation certainly depends on the class of functions 
chosen to approximate the spectrum of sample under investigation. 

It is obvious that the obtained evaluations of energy depend also on the 
precise definition of the class of events appertained to the sample. E.g. for 
the sample of penetrating events of fixed number of muons these evaluations, 
in general, depend on the size of installation (or its sensitive area) in which 
the registration of fixed number of muons is demanded to attribute the event 
to the indicated sample. 

With account of these restrictions the described method allows the re- 
liable estimation of muon mean energies by means of accompanying burst 
mean value which can be measured in experiment. 

Note that for the nontrack devices the concept of the "mean burst" can 
be substituted by e.g. "mean number of changed-state cells", or something 
like it, according to applied detecting elements. Of course this will require 
the corresponding correction but in principle differs by nothing from the 
above said. 
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Figure Captions. 

1. The calibration parameter B dependence on the exponent a value for 
the installation considered. 

2. The average energy E sensibility to the exponent a variations. Only 
systematic erreos are shown. 

3. The exponent 7 and average muon energy E dependence on the mean 
burst m of the sample. 
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Consolidated Characteristics of the Device 



\ V 

M\ 


[1] 

a 


[2] 

P 


[3] 

T, GeV 


[4] || 
K, GeV 




Consolidated Parameters f v of the Device 




0.90 ±0.05 


0.0088 + 0.0016 


35.1 + 1.7 


135.6 + 5.3 




Covariation Matrix D^ v 


[1] 
[2] 
[3] 
[4] 


+2.025 • 10~ 3 
-7.153- 10- 5 






-7.153- 10~ 5 
+2.640 • 10~ 6 
-4.863 • 10~ 4 
-3.572 • 10~ 4 




-4.863 • 10~ 4 
+2.759 
+2.227 




-3.572 • 10~ 4 
+2.227 
28.43 



Table 1: 



14 



